Spontaneous quantum emission from analog white holes in a nonlinear optical medium 
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We compute the spectrum of quantum vacuum emission from analog white holes in nonlinear 
optical media using a microscopic model fully including the effect of material dispersion. Analogies 
and differences with respect to the prediction from the analogy with quantum field theory in curved 
spacetime are discussed. Critical comparison with available experimental data is made: Even if a 
qualitatively good agreement appears to be found, further experimental investigations are required 
to rule out relevant competing effects. 
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The most celebrated example of spontaneous parti- 
cle creation from vacuum fluctuations was predicted by 
Hawking [1, 2] in the context of quantum field theory in 
curved spacetime and consists of the emission of a ther- 
mal radiation from the horizon of black holes. In the last 
decades, the extreme difficulty of detecting this emission 
from astrophysical black holes has stimulated the inves- 
tigation of analogous phenomena in condensed-matter or 
optical systems [3]. 

The first claim of observation of analog Hawking radi- 
ation in a laboratory was indeed made in [4, 5] : following 
the proposal in [6] , a strong infrared pump pulse was sent 
through a nonlinear medium (fused silica in the quoted 
experiment) and created a moving perturbation of the re- 
fractive index. As a result, the speed of optical photons 
inside (outside) the pulse results smaller (larger) than the 
pulse velocity: Seen from the reference frame comoving 
with the pulse, the leading (trailing) edge of the pulse ap- 
pears then as the analog of a black (white) hole horizon. 
Even if the interpretation of the experimental results as 
Hawking radiation is still considered as controversial by 
some authors [7, 8], alternative explanations are based 
on simplified models [9, 10]. 

A first step towards the construction of a complete 
theory of these phenomena was made in [11, 12], where 
a microscopic one-dimensional model of this system was 
developed and used to predict the emission spectrum for 
black hole configurations: In particular, the crucial role 
of the material dispersion was pointed out. The purpose 
of this Letter is to contribute to the on-going debate by 
applying the model of [12] to the analog white hole hori- 
zon at the trailing edge of the pulse from which most of 
the quantum vacuum radiation is expected to be emitted. 
As a consequence of nonlinear effects during propagation, 
this edge quickly becomes substantially steeper than the 
leading one. 

As a first step, we consider the idealized case of a large 
refractive index modulation and we highlight the con- 
sequences of the multi-branched Sellmeier dispersion of 
fused silica as compared to a simpler dispersion as the 
one of diamond. We then move to more realistic cases 



inspired by the experimental parameters [4, 5] and we 
provide predictions for the emission spectrum that are 
in a good qualitative agreement with the experimental 
results. We conclude the Letter with a brief discussion 
of possible competing effects and of experimental signa- 
tures that would provide an unambiguous proof of the 
quantum nature of the emission. 

Theoretical framework — Following [12], matter-light 
interaction in our system is described in terms of a gen- 
eralized Hopficld model [13] where the electromagnetic 
field is coupled to several matter polarization fields de- 
scribed as uniformly distributed charged harmonic oscil- 
lators with elastic constants and frequencies (see 
the supplemental material). With suitable choices of the 
parameters, this model is able to reproduce the correct 
Sellmeier dispersion [14] 
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of transparent dielectrics. For fused sil- 

ica, n = 3 resonances are needed with 
/3i,2. 3 = 0.07142, 0.03246, 0.05540 and aOi )2 ,3 = 
0.1253, 10.67, 18.13 eV [15]. 

As it is sketched in the top panel of Fig. 1, the increase 
of the refractive index caused by the strong pump pulse 
is described as a discontinuous change of the local values 
of and Qj. For the sake of simplicity, we focus on 
the steep trailing edge (located at X = vT) of a pulse 
propagating rightwards at speed v > 0: This is expected 
to be an accurate approximation since light scattering 
and quantum vacuum emission at the much smoother 
leading edge is generally negligible. Outside the pulse 
(i.e. for X < vT), and tti coincide with the tabulated 
values of the material, while inside (i.e. for X > vt) 

il' = (l + e)- 1 / 2 ^,. 



# = (l + e)A, 



(2) 



In the actual calculations, the value of e is chosen in a way 
to give the desired value of the refractive index change 
for photons in the optical window. 

Assuming that for a sufficiently long time after a ini- 
tial steepening stage the pulse shape does not appreciably 
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FIG. 1. Top panel: Sketch of the refractive index profile for 
a right-going pulse generating an analog white hole horizon 
at x = 0. Central row: single-resonance dispersion relation 
seen from the comoving reference frame in the external (left) 
and internal (right) regions. Solid (dashed) curves represent 
positive- (negative-) norm branches. Parameters flo = ^2, 
Po — 02, e = 0.3, and v = 0.66c. Bottom left panel: co- 
moving frame emission spectrum for the two positive-norm 
outgoing modes: Solid (dashed) line corresponds to the open 
(filled) dot in the dispersion. Bottom right panel: total emis- 
sion spectrum in the laboratory frame. The arrow indicates 
the wavelength at which the phase velocity equals the pulse 
velocity, i.e. to = 0. 



change during propagation, the system can be described 
in the frame comoving with the pulse at velocity v by 
a time independent Hamiltonian (throughout the Let- 
ter, quantities measured in the laboratory vs. comoving 
frame are denoted by upper vs. lower case symbols). The 
held equations can then be solved in Fourier space with 
respect to time: within each homogeneous region, the 
dispersion is given by the Lorentz transform of Eq. (1) 
and matching conditions have to be imposed at the dis- 
continuity. This procedure provides the scattering matrix 
relating incoming (i.e. directed towards the discontinuity, 
as seen from the comoving frame) and outgoing modes. 
As this matrix mixes modes with positive and negative 
norm at a given comoving frequency oj, destruction and 
creation operators are correspondingly mixed. As usual, 
this leads to a finite intensity of spontaneous quantum 
vacuum emission even for an initial vacuum state with 
no incoming particles. Once the emission spectrum in 
the comoving frame is known for all modes, the one in 
the laboratory frame is obtained by a Lorentz transfor- 
mation of each flux and then by summing over all the 
outgoing positive- and negative-norm modes contribut- 
ing to a given laboratory frame frequency fl. 



Simplified dispersion relation — As a first application, 
we study the case of a non-polar dielectric such as dia- 
mond, where the only resonances in the Sellmeier disper- 
sion (1) are in the UV and the dispersion is regular from 
Q = 0, where it is almost linear, all the way through 
the IR and the visible range. Without loss of generality, 
the calculations have been performed with a single UV 
resonance at H = O2 and /3 = fa. 

In particular, the pulse speed v is chosen in a way to 
have an analog white hole horizon. The dispersion rela- 
tion of the lower polariton branch seen in the comoving 
frame is plotted in the two panels on the central row 
of Fig. 1 for respectively the external (left) and inter- 
nal (right) region; the upper polariton branch lies in the 
UV far outside the field of view and plays no role in the 
physics under examination. Solid (dashed) lines indicate 
positive (negative) norm modes. The white hole nature 
of this configuration is visible in that no light can prop- 
agate from the horizon into the internal region. At all 
frequencies a negative norm mode incident on the hori- 
zon is present, giving rise to quantum vacuum radiation. 
The usual Hawking process is however active only below 
the maximum frequency w max . 

The emission spectra in the comoving frame are shown 
in the left lower panel of Fig. 1. The solid line corre- 
sponds to the Hawking-likc emission in the mode indi- 
cated as an open dot on the dispersion curve: as expected 
on the basis of the analogy with quantum field theory 
on curved space-times, its spectrum displays the usual 
thermal-like 1/ui divergence at low frequencies [3, 16, 17] 
but vanishes at uj > w roax - On the other hand, the emis- 
sion on the second outgoing mode (closed dot on the 
dispersion) is defined for all values of ui but its inten- 
sity (dashed line) is much weaker. By energy conserva- 
tion, the emission on the outgoing negative norm mode 
(not shown) is given by the sum of the two emissions on 
the positive norm modes. The right lower panel shows 
the total emission spectrum as observed in the labora- 
tory frame. The thermal peak of the Hawking radiation 
around u> ~ is Doppler-shiftcd to a finite Q = \"fvko\ 
with fco defined by the intersection of the dispersion re- 
lation with the lo = line. In physical terms, this corre- 
sponds to the wavelength Ao at which the phase velocity 
of light equals the pulse velocity (indicated by the arrow 
in the figure). 

Full Sellmeier — We now proceed to study the com- 
plete Sellmeier dispersion (1) where the features of quan- 
tum vacuum emission are expected to be dramatically 
modified by dispersion. We start by considering the case 
of a large refractive index jump Sn = 0.12. Even though 
this value is too large to be realistically reached in fused 
silica, still it helps understanding the basic physics. 

The three resonances at Qi split the dispersion rela- 
tion into four polariton branches. The highest one at 
O > Q3 is never involved in the physics under investiga- 
tion. In what follows, we shall refer to the other three 
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FIG. 2. Upper panels: Optical branch (f2i < SI < S^) of 
the dispersion relation of fused silica (1) as seen from the 
comoving frame. The left (right) panels refer to the exter- 
nal (internal) region of an analog white hole. Solid (dashed) 
curves represent positive- (negative-) norm branches. Param- 
eters: v — 0.66c and e = 0.3, giving a strong refractive index 
modulation Sn ~ 0.12 in the visible range. Lower left panel: 
emission spectrum in the comoving frame. Solid (dashed) line 
refer to the emission into the positive norm modes indicated 
by open (filled) dots in the dispersion shown in the upper left 
panel. Bottom right panel: total emission spectrum in the 
laboratory frame. The arrow indicates the wavelength Ao at 
which the phase velocity equals the pulse velocity, i.e. u = 0. 



as lower (SI < Oi), optical (fii < SI < SI2), and up- 
per (SI2 < H < O3) branches. When the pump pulse 
is on the optical branch [11, 12], the scattering ampli- 
tudes for the lower and upper branches are negligible. 
Therefore, in Fig. 2 we focus our attention on the optical 
branch only: the main difference as compared to Fig. 1 



is the opening of a gap for oj < 



in the external re- 



gion (left panel). This apparently minor difference has 
a dramatic effect as it removes two of the four solutions 
which arc involved in the Hawking emission process at 
low lj. As a result the emission spectra plotted in the 
lower left panel of Fig. 2 no longer show the typical 1/ui 
thermal-like divergence of Hawking emission but the to- 
tal emission goes to a constant value as lj tends to 0. A 
white hole horizon is present only in the frequency range 
(w m i n , w max ): A remaining signature of the Hawking ef- 
fect is the rapid growth of the emission for decreasing lj 
within this window. With respect to the standard Hawk- 
ing case of Fig. 1, this spectral distortion is responsible 
for the emission spectrum in the laboratory frame being 
peaked at a wavelength slightly shifted from the one Ao 
where the pulse velocity would equal the phase velocity 
(indicated with an arrow). 

Weak refractive index perturbations — As a final ap- 
plication, we investigate the case of a weak refractive 
index modulation Sn w 0.001 recently considered in the 
experiments [4, 5]. With such a small value of Sn, having 
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FIG. 3. Laboratory frame emission spectra in the case of a 
weak refractive index modulation obtained with a 1055 nm 
laser pulse. Upper panels: Gaussian pulses propagating at 
v = vq ~ 0.684c (left) and 5 = 7° Bessel pulses prop- 
agating at vg = 0.689c (right). The different curves re- 
fer to different values of the refractive index modulation 
Sn « 1.6,1.3,0.88,0.33 x 10~ 3 (solid, dashed, dotted, dot- 
dashed line). Lower panels: fixed value of Sn — 1.6 x 10~ 3 
and different pulse velocities. Left: v — 0.684c, 0.685c, 0.686c 
(solid, dashed, dotted line; S — Gaussian and S — 3,5 Bessel 
pulses). Right: v = 0.688c, 0.689c, 0.690c (solid, dashed, dot- 
ted line; i.e. S = 6,7,8 Bessel pulses). In each panel, the 
arrows indicate the wavelength Ao at which the phase veloc- 
ity equals the pulse velocity. The vertical lines indicate the 
wavelengths at which resonant radiation can occur [18]. 
The line style coding follows the one of the spectra. In the 
upper panels, Ao and are the same for all curves. In the 
left upper panel, Ao = 497 nm falls outside the field of view. 



an analog horizon requires a very fine tuning of the pulse 
velocity. As the range (w m i n ,w max ) is very narrow [12], 
we expect that the resulting spectra are not significantly 
affected by the presence or not of the horizon. 

The experimentally accessible emission spectra seen 
from the laboratory frame are plotted in Fig. 3 for differ- 
ent values of the pulse speed v and the refractive index 
modulation Sn. In experiments, the former parameter 
can be tuned cither by using pump lasers of different 
wavelengths and therefore different group velocities vq or 
by focusing the pump beam with a conic lens so to pro- 
duce a so-called Bessel pulse [19-21]. This propagates in 
fact with a higher effective velocity v s = Vq / cos S > vq 
that can be geometrically controlled acting on the Bessel 
angle S. On the other hand, the amplitude of the refrac- 
tive index modulation Sn is controlled through the inten- 
sity of the pump pulse via Sn — n-i I, n-i being the optical 
Kerr nonlinearity of the material. In the case of fused sil- 
ica, n 2 ~3x 10~ 16 cm 2 /W, which means that an inten- 
sity / w 3 x 10 12 W/cm 2 is required to get Sn m 0.001 [4]. 
Remind that the damage threshold of the material im- 
poses an upper limit to the Sn value that can be achieved. 

In Fig. 3, we consider a pump pulse of wavelength 
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A = 1055 nm, for which vq = 0.684 c. The upper 
left panel refers to the case of a Gaussian pulse with 
v = vq, while the upper right panel is for a S = 7° Besscl 
pulse with vg « 0.689 c. In each panel, four different 
values of the refractive index perturbation are consid- 
ered, Sn « 1.6, 1.3, 0.88, 0.33 x 10~ 3 (solid, dashed, dot- 
ted, dot-dashed line), corresponding to input energies of 
E = 1280, 1080, 720, 270 AiJ in the experiment [4]. In 
the Gaussian pulse case, an analog white hole horizon is 
present for all considered values of E and the dispersion 
has the same shape as shown in the upper panels of Fig. 2. 
On the other hand, for the Bcssel pulses the speed is so 
large that the dispersion has on both sides the qualita- 
tive shape of the upper right panel of Fig. 2. As a result, 
there is no (w m i n , w max ) range where an analog white hole 
horizon is present. In spite of this, negative norm modes 
are still present, leading [22] to a finite quantum vacuum 
emission. In order to perform a critical comparison to 
experimental data, some scattering process must [11] be 
invoked to justify the emission at 90° in the experiment. 
Assuming this scattering to be wavelength-independent, 
the qualitative agreement of our predictions with the re- 
sults of Figs. 4 and 3 of Ref. [4] is quite good for both 
Gaussian and Bessel pulses. 

The main difference with the experimental data con- 
sists of an appreciable shift in the position of the peak 
emission, especially in the Bessel case. A possible expla- 
nation of this feature is offered by the strong sensitivity 
of the spectra on the pulse velocity. This strong depen- 
dence is illustrated in the lower panels of Fig. 3, where 
the spectrum is plotted for various velocities ranging from 
v = 0.684c to 0.690c at a fixed value of 5n = 1.6 x 10" 3 . 
Even though a horizon is present only for the curve at 
lowest v, the transition from a horizon regime to a hori- 
zonless one is perfectly smooth. In addition to the shift 
of the emission peak for increasing v, the only qualita- 
tive change is a significant broadening of the emission. 
The nominal experimental value of the Bessel angle is 
(5 = 7° and agreement with the theoretical curve requires 
a slight shift of less than 1° (which corresponds to a 10 
relative change in o). A further, unexplained difference 
concerns the emission width, that appears to be some- 
how narrower in the experiment in the Bessel case; the 
broader experimental spectrum in the Gaussian case may 
be a result of the variation of the pulse velocity during 
propagation. 

Finally, it is worth noting that the position of the emis- 
sion peak in the theoretical spectra is significantly shifted 
from the arrows indicating the wavelength values Ao for 
which the pulse velocity equal the phase velocity. Two 
effects can be invoked to justify this shift. First, for small 
values of 5n, the emission spectra in the comoving frame 
(not shown) are generally not peaked at u = [23]. Sec- 
ond, strong additional distortions are introduced by the 
Dopplcr transformation to the laboratory frame. 



Discussion — The main conclusion of this Letter is 
the recognition of the fundamental role played by the 
material dispersion in fixing the spectral properties of 
spontaneous vacuum emission from analog white holes in 
nonlinear optical systems. Even if our one-dimensional 
model cannot provide a reliable estimate of the absolute 
number of photons produced in an experiment, still it is 
expected to accurately describe their spectral distribu- 
tion. As expected from quantum field theory in curved 
spacetime the Hawking emission seen from the comoving 
frame has a thermal-like shape in non-polar dielectrics 
when the refractive index modulation is sufficiently large. 
In other materials, like the fused silica used in the ex- 
periments [4, 5], the spectrum is significantly modified 
because the white hole horizon can only be active above 
a lower frequency threshold. Doppler transformation to 
the laboratory frame introduces further deformations to 
the emission spectrum. 

Even though our predictions for the emission spectra 
are found to qualitatively reproduce the experimental ob- 
servations, further investigations are required to exclude 
other effects, which may produce similar spectra possi- 
bly with stronger intensity. In particular, the resonant 
radiation effect observed in [18] leads to a coherent emis- 
sion at any wavelength A, sharing with the laser pump 
the same comoving frequency w/, irrespectively of the 
norm of the corresponding mode. For the experimental 
configurations of [4, 5], the locations of A; (indicated 
by the straight vertical lines in Fig. 3) almost coincide 
with the peaks of the predicted quantum vacuum emis- 
sion. Therefore, given the huge intensity of the laser 
pump as compared to the quantum vacuum emission, a 
measurement of the emission spectrum alone can hardly 
exclude this latter effect from being the most likely origin 
of the emission observed in [4, 5]. Of course, detection 
of non-classical correlations between pairs of modes with 
opposite comoving frequencies would constitute an un- 
ambiguous signature of the quantum vacuum nature of 
the observed emission. 
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